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Abstract 

We use linear response analysis and the fluctuation-dissipation theorem to derive the energy 
loss of a heavy quark in the SU(2) classical Coulomb plasma in terms of the / = 1 monopole and 
non-static structure factor. The result is valid for all Coulomb couplings T = V/K, the ratio of 
the mean potential to kinetic energy. We use the Liouville equation in the collisionless limit to 
assess the SU(2) non-static structure factor. We find the energy loss to be strongly dependent on 
r. In the liquid phase with r ~ 4, the energy loss is mostly metallic and soundless with neither 
a Cerenkov nor a Mach cone. Our analytical results compare favorably with the SU(2) molecular 
dynamics simulations at large momentum and for heavy quark masses. 
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I. INTRODUCTION 



Parton energy loss at RHIC is widely viewed as a way to probe the properties of the 
medium created during the first few fm/c of the collision. The medium is suspected to be a 

a— 

There have been a number of calculations involving parton collisional [2J, [3J, |4;_5( and 
radiative energy loss at RHIC with the chief consequence of jet quenching [8(. The 

measured jet quenching at RHIC exceeds most current theoretical predictions, most of which 
are based on a weakly coupled quark-gluon plasma (wQGP). 

The QCD matter probed numerically using lattice simulations and at RHIC using heavy 
ion collisions, is likely to be dominated by temperatures in the few T c range making it de facto 
non-perturbative. Non-perturbative methods are therefore welcome for analyzing the QCD 
matter conditions in this temperature range. An example being the holographic method as 
a tool for jet quenching analysis [9|, 1 1 CM ] - 



In this letter, we follow the approach suggested in [ll|, [12, Q, Q to model the strongly 
coupled quark and gluon plasma, by classical colored constituents interacting via strong 
Coulomb interactions. This model has been initially analyzed using Molecular Dynamics 
(MD) simulations mostly for the SU(2) version with species of constituents (gluons). The 
MD results reveal a strongly coupled liquid at T « 4 the ratio of the mean kinetic to 
Coulomb energy (modulo statistical fluctuations). The fractional energy loss is also found 
to be considerably larger than most leading order QCD estimates. 

Here, we will provide the analytical framework to analyze the MD simulation results for 
partonic energy loss in the cQGP. In section 2, we outline a formal derivation of the energy 
loss in the cQGP for arbitrary values of T. In section 3, we use linear response theory and 
the fluctuation-dissipation theorem to tie the energy loss to the non-static colored structure 
factor. In section 4, we derive explicitly the non-static structure factor using the Liouville 
equation. Some useful aspects of the plasmon excitations in the cQGP are discussed in 
section 5. In section 6, we analyze the energy loss for both charm and bottom for T =2,3 
and 4 in the liquid phase and compare them to the recent SU(2) MD simulations [15]. In 
section 7, we discuss the relevance of our results to RHIC and holographic QCD. 
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II. ENERGY LOSS 



Consider an SU(2) co 
coupled colored plasma 



ored particle of charge q a travelling with velocity v in the strongly 



The equation of motion of this extra particle in phase space 



follows from the Poisson bracket 



^ = -{H, Pl } = q a -E? n (11.1) 
with the longitudinal colored electric field 



E in = -V£ T^L = ~ V ^(t, r) (11.2) 



We note that in 



111 



the SU(2) plasma is considered mostly electric with massive constituents 



m(3 ~ 3. As a result the transverse electric contribution is absent in (111.21) . Also, (111.11) does 



not involve the magnetic part of the Lorentz force for the same reasons. The latter is 
irrelevant for the energy loss per travel length r = vt 



dK 

— =vq a E« n (t,r = vt) (11.3) 
ar 

even in the ultrarelativistic case since the magnetic force does not perform work. 

The induced colored Coulomb potential $>i n d follows from the total colored potential $ to t 
through 



k) = k) + k) = ^fcff (ii.4) 

The last relation defines the longitudinal dielectric constant with <& ex (u, k) = ^2nq5(u) — k- 
v), the colored potential caused by the extra particle in the probe approximation (ignoring 
back reaction). Thus 
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Using flll.2p and fill. 31) we have for the energy loss of a fast moving probe SU(2) charge 



dK 



dk 
T 



udu^s . 
kv \e L {uj + iO,k) 



> , , , , , .. , , (n.6) 

dr 7nr " 
after using the analytical property of €l{z, k) = €l{—z*, —k) which follows from the causal 
character of the longitudinal dielectric function as detailed below, fill. 61) is identical in form 
to the one derived for the Abelian one component colored Coulomb plasma in [16[], to the 
exception of the SU(2) classical Casimir q 2 in ( 111.61) . It is different in content through the 
longitudinal dielectric constant which now should be derived for a colored SU(2) Coulomb 
plasma. Our derivation is fully non- Abelian in the probe approximation. 




FIG. 1: Static structure factors for T = 2,3, 4 



Below we show that for the SU(2) colored Coulomb plasma at strong Coulomb coupling, 
f)II.6j) reads 



dK 
dr 



dk k 2 D S 01 (k) 
F l-Soi(fc) 



+kv 

-kv ye^u + tO,^ 



(11.7) 



with k 2 D the SU(2) Debye wave number squared. Here 



eiO, k) 



nc D1 (k)W(z/uj T ) 



(II.8) 
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FIG. 2: Soi(g): molecular dynamics simulation (a) and analytic (b). See text. 



with the thermal frequency lot = Vrk and velocity vt = y/T/m and 



W{z/ut) 



2,71 



dt 



t — zJut 



-t 2 /2 



(11.9) 



The 1=1 static structure factor S 



01 



17| 



Soi(fc) 



N 



3=1 



;n.io) 



satisfies the generalized Ornstein-Zernicke equation 



Soi(fc) 



(11.11) 



1 - nc D1 (k) 

in the colored Coulomb plasma with 1 species density n = N/V. In Fig. [1] we show analytical 
results for (III. 101) around the liquid point [17| • In Fig. [2^ we show the behavior of (III. 101) 
using SU(2) molecular dynamics simulations with the dimensionless wavenumber q = k aws 
where aws is the Wigner-Seitz radius through 1/n = 4:7ia ws /3. In Fig. [2b we show the 



analytical results for the same range of T in 
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The I = 1 contribution e± plays the role of a generalized longitudinal dielectric constant 
in the SU(2) Coulomb plasma. Indeed, for weak Coulomb coupling r < 1, — ncni ~ k 2 D /k 2 



so that Sqi ~ k 2 /(k 2 + k 2 D ). The energy loss (III. 7[) reduces to (111.61) with — ► ei. At weak 



coupling e± in (III. 8 [) is the standard Vlasov dielectric function in [161 ] . The only difference is 
in the SU(2) Debye wave number. 

III. LINEAR RESPONSE 

To construe the longitudinal dielectric constant for the SU(2) Coulomb plasma we will 



ma 
in 



te use of the Liouville kinetic equations for the time dependent structure factors derived 



17| . For that we recall that in linear response, the induced color charge density p", 



ind 



V ■ E^ nd /4:TT ties with the external potential Q b ext through the retarded correlator 



Pld(t,r) = z J dt'dr' (R(r (t,r)J b (t',r>))) <S> b ext (t' , r') (mi) 
where Jg are the pertinent color charge densities. In Fourier space we have 



4tt 

*Snd = -p Ag( W ,fc)<&L* (III.2) 

with 



A%(uj,k) = -iJ e ~ Mk -r <R(J a (t,r)j|;(t',r'))> (IIL3) 
A comparison of (1III.2j) with fill. 41) yields 



^ ) - 1 ) i ' b -% A ^ < IIL4 » 

which defines the longitudinal dielectric constant. 

The retarded correlator (1III.3P is in general a quantum object, we now show how to extract 
it from the correlations in the classical and strongly coupled SU(2) colored Coulomb plasma. 
For that, we note that the colored charge density in the SU(2) phase space is 

JS(f,r) = J dQd P Q a 6f(t,r,p,Q) (III.5) 



and that the SU(2) charge-charge correlator is 



(J%(t,r)J b (t',r')) = ^5 ab J dQdQ , dpdp'Q-Q , S(t-t , ,r-r',pp , ,Q-Q') (III.6) 

where global time, space and color invariances were used thanks to the statistical averaging. 
The time dependent structure factor S = (Sf5f) was defined in [17j . Using the color 
Legendre transform of S yields 



(JS(f, r) J»(f , r')> = 6 ab J dpdp' S 1 (t -t',r- r', pp') (III.7) 

Only the I = 1 partial wave in the Legendre transform of the color part of S contributes to 
the SU(2) charge-charge correlation function. 

The fluctuation-dissipation theorem in the classical limit ties the retarded correlator A# 
in ( 1III.3jl to the Fourier transform of the classical phase space fluctuations (1III.7I) as 



9A#(o;, *0 = ^ SxCo;, k) = -5 ab ^f 9S x (z, k) (III.8) 

The last relation follows from Si(u;, k) = —2 Im Si(a>, k) between the Laplace transform and 
Fourier transform of Si with z = u + iO. 



IV. NON-STATIC STRUCTURE FACTOR 



We have shown in [18( that the 1-color partial wave of the Laplace transform of S/ obeys 
the Liouville equation 

zS^zk-pp') - j dp^zk-pp^S^zk-p^') = S ol (k;pp f ) (IV.l) 



Sq; is the I static structure factor introduced in 



17j 



S ol (k ]P p') = nf (p)5(p-p')+n 2 f (p)f (p')h l (k) (IV.2) 
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with the Maxwell-Boltzmann distribution fo(p). The structure factor hj(fe) relates to the 
standard structure factor So; (A;) by the generalized Ornstein-Zernicke equations 



- f dpdp'Soi(k;pp') = S i(k) = l + nh l (k) = (l-nc Dl (k))- 1 (IV.3) 
n J 

The self-energy kernel Ej in (1IV.1I) splits into a static and collisional contribution in each 



color partial wave I 
We note that 



18| 



Si(zk) = - f dpdp'S^pp') (IV.4) 
n J 

with I — 1 is what is needed in (IIII.8I) . For that, we solve fllV.lj) in the collisionless limit 
with 



^i(zk;pp') « —k-p5(p-p') - —k-pnf (p)c Dl (k) (IV.5) 

m m 

We recall that the SU(2) color part of the Liouville operator is a genuine 3-body force that 
only enters the collisional contribution. 18j. Inserting (1IV.5P into (lIV.ip and using (1IV.5I) 
and (1IV.3P yield in the collisionless limit 



SK^*) = ^§r / dp- f ° iP) 



et(z,k) 



z — k ■ p/m 



with 



ei(z, k) = 1 +nc Dl (k) \ dp k P ^ m f Q (p) 

z — k ■ p/m 



and 



(IV.6) 



(IV.7) 



/ 



dp ■ 



Up) 



i(l-W(^ r )) 



z — k ■ p/m 10 

If we insert flIV.6j) into (IIII.8I) and then use (IIII.4I) . we find for / = 1 



(IV.8) 



1 



Z..2 



1 



1 



e L (z,k) k 2 nc D1 (k) €i(z,k) 
Inserting (1IV.9I) into (1II.6j) yields the announced relation (111.71) . 



(IV.9) 
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V. SU(2) PLASMON 



Before analyzing the energy loss in (111.81) for heavy charged probes, it is instructive to 
discuss the zeros of the longitudinal dielectric constant ei(u, k) = in flll.8j) as they reflect 



on the longitudinal excitations in the I = 1 channel. For that, we need the behavior of W(x) 



as defined in (111.91) with x = uj/vTk for small and large ratio k/kj^. vt = >jT/m is the 
velocity of the the particles in the SU(2) heat bath. In weak coupling QCD m ~ gT, while 
in strong coupling m ~ ttT. 
In general, 



W(x) = W R (x) + z'W 7 (x) = 1 - xe~ x2/2 ip(x) + | x e~ x2/2 (V.l) 
with ift(x) = dye y2 ^ 2 the incomplete exponential function. For k <C hp or a; ^> 1, 



W(x) « -^ + ^a;e- a;2 (V.2) 
while for k 3> hp or x <C 1 



W(x) « 1 - x 2 + i^- xe ~ x2 (V.3) 



So in the long wavelength limit with k <^ k D , (111.81) expands to 



e^, k) « 1 + nCD f } ( 1 - i^e"^ ) (V.4) 



x 2 V" 2 

For small fc, ncui(k) pa Soi(A;) ~ k 2 D /k 2 whatever the coupling in the SU(2) colored plasma. 
Thus 



ei ( W) ^l-^l-i|^j 



(V.5) 



with the plasmon frequency u p = vt kr>. So for k <C k D , the zero of (1V.5|) is 



^)-^(l-z^§e-^ 2 j (V.6) 
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The SU(2) colored Coulomb plasma supports a plasmon with frequency uj p with an expo- 
nentially small width e - uj2 / 2v T k2 both at weak and strong SU(2) Coulomb coupling V. This 



result agrees with our analytic and leading kinetic analysis in the hydrodynamical limit 18] 
. The current analysis provides the non-analytic imaginary part as well. 
The high k 3> fco limit is metallic whatever T with 



£l (x) « 1 - inc m (k) xJ 7 ^ e~ x2 l 2 (V.7) 
with a metallic conductivity x = uj/ujt = uj/vxk 

a, ( w , fc) = ^- e-«W* 2 ( V.8) 

V32^ v T k 

We note that the plasmon branch disappears at high k in ( IV. 71) as the plasma turns metallic 
i.e. a collection of free colored SU(2) particles with a classical thermal spectrum. Also the 
plasmon in (IV. 51) broadens substantially at k « &d with its real part comparable to its 
imaginary part. This point causes the plasmon contribution to drop from the energy loss in 
the colored SU(2) Coulomb plasma as we show below. 

VI. CHARM AND BOTTOM LOSS 

Inserting (111.81) into (111.71) and using the explicit form (IV. II) yields the energy loss in the 
SU(2) Coulomb plasma 



dK g 2 C F ool [ kma * „ 1 
' ' dk - 



2 



dr 4tt v 2 J k 

j pv/vt 



v/vt 



/ dxe x /2 \ [ {I - nc Dl {k)) e x 12 / x + nc D1 {k) ^{x)\ +nn 2 c 2 m (k)/2 



(VIA) 



For an SU(2) probe charge after the substitution q 2 — > g 2 Cp/4:7i' with Cp the SU(2) Casimir. 
We note that ( 1VI.1I) is cutoff in the infrared by the Debye wave number since Sqi(&) ~ k 2 jk 2 D . 
So the main contribution to the energy loss in (1VI.1I) stems from the region k > ko for which 
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FIG. 3: Surface plot of —v 2 dK/drdxdq for charm and bottom. See text. 
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FIG. 4: —v 2 dK/drdxdq versus v/vt for charm and bottom quark for fixed q. See text. 

the SU(2) plasmon is too broad to contribute as we noted earlier. Most of the loss stems 
from the metallic part of the SU(2) plasma which is the analogue as rescattering against the 
free thermal spectrum explicit in flV".8[) . 

In Fig. [3] we display the integrand in ( IVI.ll) versus the jet velocity v/vt and the dimen- 



sionless momentum q = kay/s- This is a weighted plot of the longitudinal spectral function 
along the jet velocity. The two wings at small q are the two plasmons peaks, which progres- 
sively turns into the thermal distribution at larger q. In Fig. H]we show the same integrand 
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for fixed q versus the jet velocity v/vt normalized to the thermal velocity Vy. We note again 
the 2 plasmon poles around v ~ vt at small q. The vanishing of the termal distribution at 
q = follows from the extra x 2 weight arising from the denominator of (IVI.lj) for C£>i(k) ~ 
at large k. 

Since the loss is colored with only I = 1 contributing and is metallic with only k > ko con- 



tributing, we do not see colored Cherenkov radiation stemming from plasmon emission [20|, 
nor the ubiquitous Mach cone stemming from coupling to the sound mode J2I]. While the 
sound mode contributes to S in (1III.6j) it drops in the statistical averaging as only I = 1 or 



plasmon channel contributes. The energy loss in the classical colored SU(2) Coulomb plasma 
is mostly metallic with k > kp and soundless due to the color quantum numbers of the fast 
moving probe charge. 

A qualitative estimate for the energy loss follows by using Soi(fc) ~ k 2 /(k 2 + k^) and 
saturating the integrand by k > kp, 



iK -9^l I 12 r*j e -*A (VU) 




A- 



D 



The upper divergence is manifest in ( 1VI.1I) at k 3> ko since Soi(A;) ~ 1 and Cm(k) ~ k 2 D /k 2 
through the generalized Ornstein-Zernicke equation for all Coulomb couplings. The upper 
cutoff k max ~ 27 mv which is set by the maximum momentum transfer to the thermal particle 
of mass m in the rest frame of the probe particle M ^> m. Typically M is charm and bottom, 
while m ~ gT in weak coupling and m pa ttT in strong coupling for a QCD plasma near the 
critical point. For the former v/vt ~ v^fg (weak coupling) while for the latter v/vt ~ v^pn 
(strong coupling). For v/vt 1 ( j VI. 211 reduces further to 



dK g 2 C F u 2 ^2 7 m v 



For the SU(2) colored Coulomb plasma. Aside from the Casimirs, this result is analogous to 
the energy loss in the classical and Abelian Coulomb plasma HQ. 

To assess the energy loss for varying Coulomb coupling T = (g 2 C 2 /4:7i)((3/aws), we wm 
rewrite the energy loss ( 1VI.1I1 as 
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FIG. 5: Energy loss for charm (left) and bottom (right) in the cQGP: T = 2,3,4. 



-l 



x 



^['^ dxe**' 2 (^l-nc m (q))e x2 /yx + nc m (q)i;(x)^ + 7rn 2 c 2 m (q)/2 j 



(VL4) 

where g = /ca^ys and a^s is the Wigner-Seitz radius. The units for the energy loss per 
length in (TVT.4I) follows from T jaws- F° r SU(2), Cp = 3/4 for a heavy quark, and C 2 = 2 for 



thermal constituent gluons of mass m ~ nT. aws = (3/47m) 1//3 = (J^ Q 2 ^ 4x3 ) 1//3 = 0.6883/3 
for a density dominated by black-body (gluon) radiation n = 0.244(iV ( 2 — l)/(3 3 = 0.244 x 
3//3 3 . 

In Fig. [5] we show the dimensionless energy loss following from ( jVI.41) for charm and 
bottom as a function of the probe momentum jMv, for different T = 2, 3, 4 around the SU (2) 
liquid point. The numerics have been carried using the analytic structure factor of Fig. [TJ 
The energy loss is normalized to the total kinetic energy in length L, E /L = (7 — 1)M/L. 
Since the quark velocity is maintained constant, the energy loss is seen to exceed 1 for T = 4. 
The loss is very sensitive to the Coulomb coupling T in the liquid phase. 

In Fig. [6] we show the energy loss on a logarithmic momentum scale for both charm and 
bottom. The upper curve (black) is the total loss from (jVI.41) . while the lower curve (red) is 
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FIG. 6: Logarithmic energy loss for charm (left) and bottom (right) in absolute units. See text. 

just the metallic loss following from (1VI.3H . The difference is a measure of the energy loss due 
to collisions with the low momentum part of the excitational spectrum of the SU(2) plasma 
which is plasmon dominated. These are the wings shown in Fig. HI Charm and bottom 
jets with low momenta say p ~ 3 GeV experience energy loss through broad plasmons. The 
energy loss for jets with p larger than 10 GeV is mostly linear and therefore metallic. 

In Fig. [7] we compare our analytical results for the energy loss (red curve) to recent SU(2) 
numerical simulations (black curve) using the same model [ijj]. We note that the numerical 
simulations in [ijj] are quoted for the mean-potential to kinetic energy ratio V/K ~ 3 which 
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happens to fluctuate by about ±1 inside the simulation box. To make the comparison 
meaningful, it is better to use in the notation of [15] T = 1/Tay/s ~ 4.8 for n = 1/A 3 and 
A = 1/3T with A the minimum of the potential in the same notations. With this in mind, 
our analytical results at T = 4 compare favorably with the molecular dynamics simulations 
at large momenta and for heavier quark masses (say bottom). Most of the discrepancy with 
the simulations is at low momentum where the effects of the hard core in [15j are the largest. 
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FIG. 7: Energy loss: (red) analytical versus (black) SU(2) molecular dynamics 
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VII. CONCLUSIONS 



We have analyzed the energy loss of fast moving charm and bottom quarks in an SU(2) 
color Coulomb plasma for a broad range of the Coulomb coupling. The Coulomb character of 
the underlying interaction retained classically make the energy loss entirely described by the 
longitudinal part of the dielectric function. We have used linear response theory to derive an 
explicit expression for the imaginary part of the dielectric function in terms of the Laplace 
transform of the time-dependent structure factor in the SU(2) Coulomb plasma. 

We have shown that the probe initial color and statistical averaging causes the longitudinal 
dielectric function to select the / = 1 color channel of the time-dependent structure factor 
which is the plasmon channel. The sound channel dominates the low momentum of the I = 
color channel, and decouples from the longitudinal part of the dielectric function. While the 
SU(2) plasmon survives at strong coupling, its width for k > ko is substantial and therefore 
causes it to thermally decay. 

The energy loss of fast moving charm and bottom quarks is mostly due to the metallic 
aspect of the SU(2) colored Coulomb plasma which is dominated by thermal particles. There 
is no colored Cerenkov cone as the plasmon is dwarfed in the metallic limit, nor a colorless 
Mach cone as the sound decouples due to the probe initial colors. The energy loss is soundless. 
Our results are of course only classical. They apply for a broad range of V near the liquid 
point. The comparison to the MD simulations show that our energy loss is about comparable 
at higher momenta and for heavier quarks where the effects of the numerical hard core is 
small. As initially reported in [15], the energy loss is sizable. 

Strong coupling assessment of jet energy loss in gauge theories have been carried out in 
the context of holographic QCD 22]. The fact that a Mach cone was reported in these 
calculations [23[], maybe due to the fact that the probe jet is actually colorless. Indeed, most 
of the holographic jets are inserted with an external hand that maintains a constant velocity 
and perhaps even balance the color charge. Clearly colorless (mesonic) jets of the type QQ 
do couple to the sound channel in our case through the S o(&) structure factor [171 ]. and are 
expected to be trailed by a Mach cone. 

Finally, to carry our analysis of charm and bottom at RHIC and perhaps even LHC, 
require an assessment of the heavy quark composition in the prompt phase of the heavy 
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ion collision which we have not carried out. Also, we need to address more carefully the 
correspondence between our classical SU(2) QGP and the quantum SU(3) QGP. These issues 
will be addressed next. 
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